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Abstract 



We calculated response functions of the deuteron for charge exchange processes, 
' including the final state interaction between two protons. Using them we evaluated 



the double differential cross section and polarization observables of d(p, n)2p by means 
of plane wave impulse approximation with an optimal factorization. Calculation well 
reproduced the shape of the energy spectra of the cross section, though somewhat 



overestimated the magnitude. It also reproduced the spin observables well. 
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1 Introduction 

Intermediate energy nucleon-nucleus scattering with large momentum transfer (|q| >, 1 fm _1 ) 
shows a broad bump in the energy spectrum. This is called a quasi-elastic scattering because 
the scattering from an individual particle in the nucleus is expected to play a dominant role 
in the process. The peak is observed near the nucleon recoil energy q 2 /2Mn(Mn being the 
nucleon mass) and its broad width is thought to reflect the Fermi motion. 

For an external field which induces charge exchange processes, the nuclear response is 
separated into those for isovector spin-scalar r, isovector spin- longitudinal r(cr • q), and 
isovector spin-transverse t(<t x q) modes. We write the response functions for these modes 
as Rs, Rl, and R? respectively. It is a very interesting subject of recent nuclear physics to 
extract them from the quasi-elastic scatterings. In the early 80's Alberico et al.[[I[ predicted 
that the precursor phenomena of the pion condensation are prominent in this region. Their 
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calculation showed that in the high q transfer region Rl is enhanced and Rt is quenched 
relative to those in the free Fermi gas. Therefore the ratio Rl/Rt should be much larger 
than 1. 

Since the polarimeter technique has improved greatly these days, complete polarization 
transfer experiments have become feasible. Now the observables like the depolarization 
tensor Dij as well as the analyzing power A y are available. The prediction of Alberico et al. 
inspired experimentalists to extract Rl and Rt from such polarization observables. 

The transverse response function Rt has been measured in a wide range of (a>, q) by 
electron scatterings Q for a long time(u; being the energy transfer). However we need to use 
the hadron probe instead to excite the spin- longitudinal mode. The first experiment which 
extracted the response in this mode was performed by Carey et al.@ at Los Alamos meson 
physics facility(LAMPF). They carried out the inclusive (p,p') scattering experiment using 
500 MeV proton. The targets were 2 H, Ca, Pb. However they got the result of Rl/Rt £ 1) 
which contradicted the prediction. 

The (p,p') reaction does not distinguish the isoscalar response from the isovector one, 
while the (p,n) reaction makes a good probe of Rl because of its purely isovector nature. 
Lately Chen et al. ||] performed polarization transfer measurements for the quasi-elastic (p,n) 
reaction at LAMPF with the incident energy 495 MeV at the laboratory angle $i a b = 18°. 
The targets were 2 H, 12 C, 40 Ca. Further experimental studies are now under way for 
various angles and bombarding energies at LAMPF and the Research Center for Nuclear 
Physics, Osaka University(RCNP) §. 

In the analysis of these experiments the data of the deuteron played the essential role, 
the role of reference. When Carey et al. derived the ratio Rl/Rt, they assumed R d L /Rr — 
l(d referring to the deuteron) regarding the deuteron as an assembly of a free proton and 
neutron. However there is a possibility that the correlations among constituent nucleons, 
such as the well-known tensor correlation, are important for determining the ratio. Very 
recently Pandharipande et al.|j) calculated the response functions R d L and R T including such 
correlations in the initial deuteron and final 2p states by using Urbana-Argonne interaction. 
They found that the ratio Rl/R t is close to 1 at the peak region but deviate rather much 
from 1 at high and low ui. In this paper we have calculated them together with R d s making 
use of the Reid soft core potential JF). We have obtained the same result as theirs on the 
ratio(see sect. 4). 

We also need to know the reaction mechanism for extracting the response functions from 
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the experimental data. We have considered d(p,n)2p reaction, because the deuteron gives a 
good test case for that purpose. Since it is a weakly bound two nucleon system, we expect 
simple mechanism for such intermediate energy reacitons. Furthermore reliable calculations 
are feasible in this nuclear system. 

Early theoretical studies of this reaction were done by Phillips ||, Dass and Queen 
0. They calculated D nn at the incident energy of 30 ~ 160 MeV with inclusion of only 



S-wave correlation of the final 2p. Ramavataram and Ho-Kim |10] treated the reaction 
at the incident energy of 25 ~ 100 MeV with different sets of NN phase shift. As the 
polarization measurements improved in the intermediate energies, theoretical studies also 
made progress. In this energy region Bugg and Wilkin investigated the tensor polarization 
of the p(d, 2p)n reaction where the outgoing two protons has very low excitation energy [p| . 
Carbonell, Barbaro and Wilkin also considered the same reaction in a more advanced scheme 
12]| . Recently Deloff and Siemiarczuk |13] considered relatively high energy scattering by 



the impulse approximation. They treated pd — > pnp reaction at 1 GeV including the NA 
channel in AW scattering. 

Here we have taken the plane wave impulse approximation(PWIA), and followed the 
optimal factorization formalism described by Ichimura and Kawahigashi ]14| and the method 



of partial wave development of Carbonell, Barbaro and Wilkin ||12|| . We have calculated the 
unpolarized double differential cross section and the polarization quantities of the reaction 
d(p,n)2p performed by Chen et al.0 with the response functions mentioned above. 

In sect. 2 we present the formalism to calculate the i-matrix of the reaction in PWIA 
with the optimal factorization. In sect. 3 we define the response functions and relate them 
to polarization observables. In sect. 4 we give the results and the discussions. 



2 Formalism 

2.1 ^-matrix of d(p, n)2p 

In the impulse approximation we take the transition operator of this reaction T as the sum 
of the AW transition operator toi- Here we refer to as the incident particle and i as the i-th 
nucleon in the target. The deuteron is a weak binding system and its density is comparatively 
low. In addition the projectile has kinetic energy of the order of several hundred MeV. So 
one may expect that the PWIA is reasonable for this target, the confirmation of which is 
one of our main purposes. 
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The t-matrix element of the reaction can be written in PWIA as 

(f\f\i) = J dk 1 dk' 1 6®(k f + k\ -k,- k ± ) J dk 2 dk' 2 S i3 \k 2 - k' 2 ) 

x (x^f', n \(ff\k'ik , 2)(k f k'i\io 1 \k i k 1 )(k 1 k 2 \fi)\x^;p) 
+ (1 * — >2) 

= 2 J dr x dr\\ dr 2 dr' 2 5®(r 2 -r' 2 ) J * ^(rVrO^r^n 



x 



(27T)* 

(<p f \r' ir ' 2 )(k f p - kfltoxlkt p - k^ir^l^lx^py^ (2.1) 



where tfi and cpj are the initial and final states of the target, \Xm'iP) an d \Xfj. f i n ) are the 
spin and isospin states of the incident proton and the outgoing neutron respectively. We 
separate the center-of-mass(c.m.) motion from (r 1 r 2 |v?i) and {r'xr' 2 \ipf) as 

( ri r 2 | Vi > = (rm )M )e iK < R , (r'^f) = (A^ Ms )e iK t R ' . (2.2) 

Here R and r are the cm. and the relative coordinates; likewise for R' and r'. The subscripts 
of \I/ mean the total intrinsic spin and its projection. The final 2p does not have a bound state 
and is specified by the aymptotic relative momentum k and its spin (S,Ms). Performing 
some integrations we obtain 

(f\f\i) = (27f) 3 5 (3 \k l + K i -k f - K f )f fi , (2.3) 

f fi = 2jdrdr' J ^(xM/^K^MjO^/p-fe/ltoil^p-feiJ^l^^lxwJp) 

Xe i{\K f +k a -p)r e -i{\K f+ k a +\q-p)r' ^A) 

with the momentum transfer q and the average momentum k a , 

q = k f -h, k a = ^±^l. (2.5) 

It is very complicated to advance the calculation, because there is p dependence in the AW 
^-matrix, and it is off-energy-shell in general. To proceed furthur we replace p in the NN 
t-matrix with a certain fixed value p which well represent the contribution to the integral. 
Equation ( |2.3| ) becomes after the integration over r' and p 

f fi = 2j rfr"<^;^.; r^KvI/^^I^Xfc^ ^ — fc^|£ i|fc, ^» — fc,> <^ | vl/^^^) | ; ^> e ~ . (2.6) 
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From now on we will evaluate Tfi in the projectile-target cm. frame. Following the 
prescription of optimal factorization [13 \T5\ we set the optimum value p as 

k n 



P 



-j-VQ, 



where rj is determined by the on-energy- shell condition 



(2.7) 



E(k f ) + E(p - k f ) = E(ki) + E(p - ki) 



(2.8) 



with E{k) = \j Mjj + k . The NN t- matrix in Eq. ( |2.6|) now becomes 

(k a + ^q, (-g - v)q - ~K\ioi\k a - -q, (- - v)q - \k) 



(2.9) 



and we call this the NN t-matrix in the ^-frame. 



2.2 Observables 

We introduce the orthogonal unit vectors in the p-d cm. system as 

q = - — -, n=— p = q x n. (2.10) 

\q\ x «/l 

The unpolarized differential cross section J and the depolarization tensor Dij are expressed 
as 

I = A - Tr| fft]. A J = Tr '^° jl , (2.11) 

where the subscripts « and j represent one of the directions q, n and p, whereas the subscript 
indicates the projectile. In the equations above Tr[OiO|] means 

Tr^Ol] = EEE /^(Xm/^I(^)mJ0iI*(i)m)Ix W ;p) 

MiM/ M SM S 

x(x / , i ;p|(^ 1) M|0||*g ) M S )lxM / ^>^-(^ / -^)) ) 

where a; = E(ki) — E(kf), and and £^ are the target energy of the initial and the final 



states respectively. The energy of the deuteron £ l is given by y Mj + k i . We approximate 



£f as yJ(2M p + k 2 /M p ) 2 + kj, because we consider only the region where the relative motion 
of 2p can be treated non-relativistically. The kinematical factor K in Eq. ( [2 . 1 1| ) is given by 



rriimf kf 1 , . 

K = - ± -i^ 1 - 2.12 
(2vr) 2 fci 6 V ; 



where rrii and mj are the relativistic reduced masses 

E{k i )E i E(k f )£f 



nii 



(2.13) 



E{ki)+£ v J E(k f )+£f 
The unpolarized double differential cross section in the laboratory frame i/ a & is connected 



to Jq through the relation 



I i ab sin 0j a6 = I sin C 



Bleszynski et al. |T^] introduced the quantities expressed by Djj as 

J), 



■[1 + D m + D ? , + D 



[ 1 I D nn Dqq -Dpp], 



n n T~ Dqq Dpp i 



[1-D nn -D qq + D 



(2.14) 



/'/'J 



(2.15) 



We will show the results in terms of them. 



2.3 NN t- matrix 

In Eq. ( |2.6|) T/j includes the NN t- matrix in the 77-frame, which we should evaluate from 
the one in the NN cm. frame (kj, —Kf\toi\Ki, —Ki). This is described in the convention of 



Kerman, McManus, and Thaler as 

(«/, -Kf\t Q1 \Ki, -Ki) = —, — -M{Ki, Kf) 

M(K h K f ) = A + B{trV> -n c )(o-(°) • n c ) + C((t (1) -h c + tr^ • n c ) 

+£(<r< 1 ) ■ q c )(<x(°) • q c ) + F(*W ■ p c )(tr(°) ■ p c ) (2.16) 

with the orthogonal unit vectors 

Q c = H r> ™ c = I — — -77, P c = Qc x n e . 2.17 

Here we used = = |k| and M(ni,Kf) is the AW scattering amplitude. The coeffi- 
cients A ~ F are determined by the AW scattering experiments. Due to charge independence 
they are written with the isospin operator as 

A = A + Ai(t<°> • r (1) ), etc. (2.18) 
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The relation between the AW t-matrices in the ?7-frame and in the NN cm. frame is 



(k a + ^q,(~-r])q-^k a \i l\k a _ } 



p.19) 



11 A 1 % 

J{k a , q, 7]) exp(-x'cr 0n ) exp(--pVi„) (n f , -K/|*oi|«», exp(-xV 0n ) exp(--pcx ln ) 



Here J(k a , q, rj) is a Moller factor 
J(k a ,q,rj) -- 



E( K ) 



--. (2.20) 

/£(fc a + lq)E((-rj - ±)q - |fc a )£(fc a - ±q)E((- V + \)q - \k a ) 

and x, x', P, and p' are spin rotation angles defined in Ref. |14| . 

Now we introduce two approximations. One is to identify the unit vector q in the ?7-frame 
with q c because the angle ip between them is extremely small for wide range of u as is shown 
in Table I. Since we determine the optimal value by Eq.( |2.7| ), the reaction plane is common 
for the 7] frame and the NN cm. frame. Then 



P~Pc- 



(2.21) 



The other approximation is to neglect the spin rotation, because the rotation angles are very 
small in the present case as shown in Table I. The NN t- matrix element in the rj- frame now 
turns to be 

(k a + \q, (-- - r])q - -k a \i 01 \k a - \q, (~ - rj)q - ]-k a ) 



2"' v 2 2 
J(k a ,q,r))(nf, -Kf\t i\ni, 



[2.22) 



Then we express the t-matrix ( |2.6| ) as 

T f i = 2C[A 1 <x M/ ;n|rW|^;p><^J r jT« e -i9- p |«f 1)M > 



+Bx( X 



ft f , 



n\T 



{0 \<r^ ■ n)\ X ^vW {S)Ms 



,2p k 



rW(o- (1) -n)e-i^ r i* 



(1)M/ 



+c 1 { x , f M^\^.u)\ x ,-p){^ 

+C 1 ( X , f -n\r^\x,-p)(^^ 
+Ex(x 



[2.23) 



n\T 



(0)^(0) . 



P)e-^ r \<w) 



{S)M S 



-(1)^(1) 



with 



47T 

EM 



(2.24) 



where we used Egs- ffTTSp , d2T|), fl2g| ) and (^22|) . 
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2.4 Wave function of the target nucleus 

We describe the wave function of the deuteron as 



\n p) — \p n) 



(rK )M ) = J \ (2-25) 



(r\$ d im ) = E Rx(r)[Y x (n r ) x xi} ( $ (2.26) 

A=0,2 

where Y/ m (f2 r ) is a spherical harmonics and Xsm s is a s P m wave function of the two nucleon 
system with the spin S and its projection Ms- 
The scattering wave (^|^(s)m s ) * s expressed as 

(r\*Ts)M s ) = (r\*? s) M s )\PPh (2-27) 
(r\*Ts)M s ) = 4EE E i l e^^^(lrnSM s \JMj)Y; m (n k ) 

V n JMj Im /'+5=even Kj ' 

x[Y v (Q r ) x xs^], (2.28) 

where 07 is a Coulomb phase shift. The sum over V is taken for I' + S =even because of the 
identical protons. We solve the Schrodinger equation for ipffi (r) with the asymptotic form 



^i% r ) — ► ^5 U/ exp(-i(kr - - r]\n2kr + a tl )) 



-Sff exp(i(£r - - r]\n2kr + <r r )) j (2.29) 

where i] is a Sommerfeld parameter. Due to time reversal invariance the S-matrix S is 
symmetric unitary matrix which can be diagonalized as 

Sfi = Y,Uy<Ui* (2.30) 

a 

with a real orthogonal matrix U. Using the wave function ip// S a (r) = U^„i/ji, 6 i„(r), we 
rewrite (r|$g } ^ s ) as 

V" JMj a ; m i'+5=even fhl 

xilmSMslJM^Y^n^Y^r) x Xs ]<£. (2.31) 
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3 Spin Response Functions 



In the calculation of the observables ( 2.11|) by Eq.( [2.23 ) we encounter the following four 
response functions defined as 



,1 (i)m)\ 2 



Rs(q,<*) = ^EE /^l(^)M S |T« e -i^|^ 

6 M SM S 

x5(u - (£ f - £ 1 )), (3.1) 

6 M SM S 

x5(u - (£ f - £*)), (3.2) 

6 M SM S J 

xS(u-(£ f -£*)), (3.3) 

|2 



'fl)A/)| 2 



6 M SM S J 

x8{u-{£ f -£ i )). (3.4) 



They represent the excitation strength of the nucleus induced by the external field with 
4-momentum transfer (u>,q). Here Rs is referred to as the spin-scalar, Rl as the spin- 
longitudinal, and Rt as the spin-transverse response function. We will see that Rt, u = Rt, p 
later. The factor | means the average over the initial spin states of the deuteron. From 
Eqs.(2.25) and ( |2.27| ) we can easily calculate the isospin part, for example, 



From now on we use the spherical tensor representation such as 

ci±i = T(&x ± icr y )/\/2, (7io = o'oo = 1, (3.5) 

<?i±i = =F(gx ± iq y )/V2, qio = q z , etc. (3.6) 

and 

<T (1) -g = E^2C> etc. (3.7) 

Sz 

We can rewrite Eqs.( |3.1|) ~ ( |3~4] ) as 



flsfow) = (3.8) 

Jfefow) = EC<k^fe ( 3 - 9 ) 



9 



Rr, n {q,u) = Y,n\M R % ( 3 - 10 ) 

s z s' z 
s z s' z 

where 

U s's' z - ^2^2^ aK \^(l)M\ a s's' z e2 \^(S)M s /\^(S)M s \ a ss z e 2 l*(l)M/ 
° AI SM S J 

x5(u-(E f -£*)). (3.12) 
Using the expressions (|2.26 ) and ( |2.31 ) we get the matrix element above as 



12EEEEE E 2Z*- L - l (-r s 

A JMj Im LM L I Mi l'+S=even a 

x (lmSM s | JMj) (LM L ss z | JM 7 ) (1MIM, 1 JMj) (AOLO |Z'0) 

x^(2J + 1)(2L + 1)(2A + 1)(25 + l)e-"'t#frg s J (k, q)Y lm (Sl k )Yl ML (Sl q ) 



I' S J 
A 1 1 
L s I 



M\\WW\\m\s\l-±8) { A 1 1 } (3.13) 



where 



(V^llVl^ forS = ° (3 14) 

and T1\ SJ (k, q) denotes the integration of the radial part, 

rg s J (fc, g) = J r 2 dr^f^j L (-qr)R x (r). (3.15) 

Now we apply Eq. (|3.13|) to Eq. (|3.12| ). The integration over the angular part of the relative 
momentum k can be analytically performed because the argument of the delta function has 
no angular dependence due to the approximation for Sf. By use of the orthogonality of 
Clebsch-Gordan coefficients, E/^f, and spherical harmonics, R S J*, results in 



= 48EEV( 2i i + 1 )( 2i 2 + i)^f ili2 (fi ! ) 

x J k 2 dk5{u-{£f -n)Ri LuS , L2 (k,q) (3.16) 



where 



^W^) = E ^ L i M l,ss z \IM i ){L 2 M L2 s's' z \IM i ) 

M Ll M L2 M 7 
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xYZ iMLi (n q )Y L2ML2 (n q ), (3.17) 
^WM = EEE E E E* L2 - Ll (-r +s ' 

S AiA 2 J i' 1 +5=even i^+5=even « 

x(2J+ i)(as + + i)(2A 2 + i)rS Jt (M)rS J (M) 

x(A 1 0Z. 1 0|Z' 1 0)(A 2 0Z. 2 0|Z^0)<||| C xi 1 )|||>{||| C riP|||> 



.11 1 1 1 1 



^ ^ I f ^2 S J 



y\V(~S-l;-s)W(~S-l;-s')\ A a 1 1 || A 2 11 }.(3.ls) 



In this form we are able to separate out f2g-dependent part. 
When s = s' = 0, we obtain 



dOO 
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£(2/ + l) 2 / k 2 dk8(u — (£ f - Z^Roifliik, q), (3.19) 
1 J 

R&ijufoq) = p^EEEE E( 2 ^ + 1 )v / (2^+^)(2^Ti T 



A1A2 J Z' 1= oddZ^=odd a 



xrS 1J (fc,g)rS 1J *(fc,g)(A 1 0/0|/;0)(A 2 0/0|/ 2 0) 

xlf(U;i/; JA 1 )W r (l/ 2 l/; JA 2 ). (3.20) 

In the case of s — s' — 1 the factor i?{ Li 1L g) is rather complicated, so we leave it as 
it is. On the other hand we can sum up over s z and s' z in Eqs.(|3.9|), ( |3.10| ), and ( |3.11| ). For 
the longitudinal mode 



^ , , J(2Li + 1)(2L 2 + 1) 

E vtjs'T^UM = ^ £ ^LiOlOl/OJCLaOlOlJO) (3.21) 



and for the two transverse modes 



^(-) m {^(-) Ll ^,L 2 (3.22) 



+ -^=(LiOL 2 0|20)W(LilL 2 l; I2)y/(2L X + 1)(2L 2 + 1) } . 



The derivation of Eq. fl3.22| ) is in the appendix. Thus the two spin-transverse response func- 



tions (|3.3| ) and (|3.4j) are identical 



Now we write down the spin observables Q2.ll ) in terms of the scattering coefficients 



A ~ F and the spin response functions. When we calculate Iq and D i: we have only those 
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response functions in Eq.( |3.1| ) ~ Eq.( |3.4| ) and the interference term does not appear . From 
Eqs.( |2.23| ) and ( |3~1"| ) ~ ( |3.4j ) the unpolarized double differential cross section becomes 



J = 48^C 2 {(|A| 2 + |C 1 | 2 )i? s +(| J B 1 | 2 + |C 1 | 2 + \F 1 \ 2 )R T + \E 1 \ 2 R L }. (3.23) 
As for the polarization observables Di defined in Eqs.( p.l5|) , 



I D = 48K( 2 ({A^Rs + IC^Rt 

I D n = ASK^^C^Rs + IB^Rt 

I D P = ASKC'IF^Rt, 

I D q = 48^C 2 |^i| 2 ^l. (3.24) 

4 Results and Discussions 

For the deuteron we used the wave function for the Reid soft core potential [7j. We solved 
the Schrodinger equation for ipiH(r) in Eq.( 2.28|) with the same potential which is supplied 



for J < 2. For J > 3 we neglected the nuclear interaction and set ipi/i(r) = Fi(kr)5i>j, where 
Fi(x) is a Coulomb wave function. We used two different NN scattering amplitudes. One is 



provided by Bugg and Wilkin ||19|| , and the other is calculated in the SAID system from the 



phase shift SM89 of Arndt et al.fl20fl 



In Fig.l we show the results of the response functions (&)Rs, (Jo)Rl and (c)Rt- The 
solid lines show the results with full correlations of the initial and the final target state. 
The dashed lines indicate the results with the final state correlations. The results when the 
deuteron consists of only S'-wave are presented with the dotted lines. The dash-dotted lines 
are the contributions of the final S- wave to Rl and Rt- All three response functions have 
large peaks at ui ~ 60 MeV and their general features are almost the same. Comparing the 
solid line with the dashed line we see the effects of the final 2p correlations. They commonly 
reduce the quasi-elastic peak and slightly shift Rt toward lower uj. As is expected, the 
S- wave interaction of 2p forms a sharp peak near the threshold in Rl and Rt- For Rt it is 
much smaller than for Rl- This difference is due to the D-state of the deuteron. As we see 
in Fig. 1(b) and (c), the threshold peak is almost the same in Rl and Rt unless this -D-state 
exists. From this we see the interference term with the initial D-state and final S-state 
makes opposite contributions to Rl and Rt- 

Figure 2 shows the ratio Rl/Rt- Except for the dash-dotted line the meaning of each 
line is the same as in Fig.l. In the dash-dotted line we replaced the final P-waves by those 
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without correlations. The behavior of the solid line is essentially the same as the one obtained 
in Ref . || . Comparing the solid and the dashed lines we find that the ratio becomes smaller 
at the lower u side and becomes close to unity at ui ~ 50 MeV when we include full final 
state interactions. At higher u side it increases gradually up to 1.4, though the ratio remains 
stable at 1.1 in the case of the uncorrelated final state. The interference term with the initial 
.D-state and final P-state contributes differently to Rl and Rt, and the P-wave correlation 
amplifies this feature. The initial P-state plays an essential role in raising the ratio near the 
threshold. 

We display the results of the unpolarized double differential cross section in the laboratory 
frame i/ Q & in Fig.3 and the polarization quantities Dq, D n , D p , D q in Fig.4. The solid lines 
and the dashed lines represent those obtained with and without the final state correlations 
respectively. They are calculated by use of the AW amplitude of Bugg and Wilkin [p 



The dotted lines show the results with the final state correlations, but calculated by the 



NN amplitude of Arndt et al.pO|. For i^, the peak position and its width is almost the 
same as the experimental result, but the theoretical magnitude is more than 10% larger. 
Comparing the result of Bugg- Wilkin's amplitude with that of the Arndt's amplitude we 
find the uncertainty caused by the ambiguity of the AW amplitude. It is about 4 % in this 
case. So this is not the only reason for the deviation from the experimental result. A small 
peak due to the S"-wave correlation of the final 2p appears near the threshold in solid and 
dotted lines. 

The polarization quantities Dq, D n , D p , and D q are in rather good agreements with the 
experimental data in the range of errors. The calculation by the uncorrelated final state 
does not differ so much from the one with the correlations. However the differences between 
the results by the two AW amplitudes are distinct in D n and D p . The behavior of P n , D p 
and D q can be understood with a help of the ratio Rl/Rt- For example P„ is written as 



D = |W| ^g/^-ri^ii 

^ n HA 19 , I/! IWn ID , t \ D 19 , l/t 15 , I 77- m , I 771 19 D ID V*' L > 





Ci 


2 Rs/ Rt 


+ 1 


Pi 


2 




2 + 




2 )R S /R T + 


m 


2 + 




l 2 + 


Pi 


2 ) + 


Pi 


2 Rl/Rt 



from Eqs.( |3.23D and ( |3.24 ). The values \A\\ 2 and |Ci| 2 are considerably smaller than |Pi| 2 , 
|Pi| 2 , and |Pi| 2 . The numerator is almost determined by |Pi| 2 , and |Pi| 2 Pi/P-r is the main 
term which changes the denominator. From the behavior of the ratio Rl/Rt shown in Fig.2, 
D n is expected to have a broad peak at w ~ 50 MeV and a sharp decrease near the threshold. 
We can derive the analogous equations for D p and D q . They show that D p is expected to 
behave similarly as D n and D q to have a broad dip at uo w 50 MeV and a sharp rise near 
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the threshold. Calculated results have these features, though they are not apparently seen 
in the experimental data. 

On the whole our calculation agrees rather well with the experiment, but it somewhat 
overestimates 1^. One reason of this discrepancy is the ambiguity of the NN amplitude. 
We have seen that this causes the uncertainty of several percent in 1^. There are also clear 
differences in D n and D p between the results by the two amplitudes. Another reason is the 
applicability of the optimal factorization. In Table II we present the values of the relevant 
quantities. The value of rj increases rapidly as u. On the contrary the kinetic energy of the 
projectile in the NN laboratory frame T)^ decreases correspondingly. The NN amplitude 
varies rather much as a function of this . We consider that the optimal factorization-like 
treatment is needed, but there may be other ways of estimating the internal motion. A naive 
treatment of the reaction is to assume that the target nucleon is at rest in the target, which 
corresponds to the case of 77 = 1/4 as is seen from Eq.( [2.7|) . At uo = 60 MeV, 77 is 0.265. So 
this picture is partially realized near the peak. We may attribute the overestimation of Ii a b 
largely to neglecting the absorption of the flux, such as a shadow effect. Virtual A excitation 
may also be a reason. 

In summary the ratio Rl/Rt is close to unity near the peak. Carey et al. 0] used the 
quantities of the deuteron for reference by taking average of all uo. Since most of the data 
are taken near the peak, their prescription of Rl/Rt = 1 is permissible, or it should be 
slightly larger if we want to fix it for convenience. We can also conclude that PWIA is a 
reasonable approximation for the reaction d(p,n)2p. It gives rather good results for D{. It 
also well reproduces the shape of Jj a &. However it overestimates the magnitude, which cannot 
be explained only by the uncertainty of the AW amplitude. 



The numerical calculation in this work was performed on the TKYVAX node in the 
Meson Science Laboratory, University of Tokyo. This work was supported by Grants-in-Aid 
for Scientific Research of the Ministry of Education (No.02640215, 05640328). 



Appendix 

From Eq.( p,17| ) Eq.( |3.22|) is calculated as follows: 



S Z S' Z 
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A-jr 

T E E Y.( LM l^ z \IM i ){L'M li \s' z \IM i ) 

6 s z s' z M L M L , Mi 

x(-) s ^i ^n)Y ls ^n){-) ML Y L _M £ (n 9 )y L , Mi ,(o 9 

/ A — I 

-s z +s' z +M L -h+I-L' 



1 An 


/ An 




/2L' + 1 


/ 2Zi + 1 V 


2l 2 + 1 V 


An V 


An 



E E EE 

s z s' z M L M L , Afj i^ia 

xF ;i _ Sz+s ,(Q n )^ 2 _ M£+Mi ,(Q g )(1010|Z 1 0)( J L0L'0|i 2 0) 

x(L -M L L'M L i\h - M L + M L/ )(LM L ls z \IMj)(IMjl -s z \llM v 

x(ls z l - s' z \h s z - s'J 

21 + 1 



= ^(1010|/0)(L0L'0|/0)W(L1L'1; Il)———^{2L + 1)(2L' + 1) 

2/ + 1 

x(-) /+m ^_ m (fi n )^ m (^) (a.i; 

where we use the relation 



(2e + 1)(2/ + l)(aa/ 7 - a|c7)W(a&cd; e/) 
= ^(aai&/3|e at + /3)(e a + P d 7 - a - (3\cy)(b/3d 7- a- 7 - a). (A. 2) 

In Eq. (|A.l| ) we can sum over m by using the identity 

2/ + 1 2/ 4- 1 

E^mW^W = ■ g) = -^-^(O)- ( A -3) 

Due to the presence of (1010|Z0), we have only to take the sum over / = 0, 2. After inserting 
the explicit values of Pi(0) we obtain Eq.( j3.22 ). 
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Table I. which is the angle between q and q c and 
the spin rotation angles Xi X' ■> P an d P' ■ 





40 


60 


80 


100 


120 


MeV 


i> 


0.072 


0.222 


0.363 


0.492 


0.605 


deg 


X 


0.936 


-0.016 


-0.911 


-1.726 


-2.440 




x' 


1.826 


2.279 


3.579 


4.353 


5.028 




p 


0.967 


-0.017 


-0.933 


-1.756 


-2.464 




p' 


1.891 


2.824 


3.697 


4.481 


5.152 





Table II. The values related to the optimal factorization: 
7], the corresponding incident energy in the NN labora- 
tory frame T^ b N and the scattering angle 9™. 





40 


60 


90 


120 


MeV 


V 

rpNN 

A lab 

oNN 

a lab 


0.090 
513 
17.7 


0.265 
493 
18.0 


0.519 
456 
18.9 


0.751 
413 
20.2 


MeV 
deg 



Fig.l. Response functions (a) Rs, (b) Rl and (c) Rt- The solid line includes full correla- 
tions. The dashed line has no correlations in the final state. The dotted line is the result of 
purely S'-wave deuteron. The dash-dotted line indicates the contribution of the final S'-wave. 

Fig. 2. Ratio Rl/Rt- Lines are the same as in Fig.l except for the dash-dotted line. The 
dash-dotted line has no correlations in the final P-waves. 

Fig. 3. Unpolarized double differential cross section in the laboratory frame Ii a b- The solid 
line includes full correlations and the amplitude of Bugg and Wilkin is used. The dashed 
line does not include final state interactions and the amplitude of Bugg and Wilkin is used. 
The dotted line includes full correlations and the amplitude of Arndt et al. is used. The 
experimental data in Ref. || are represented by the dots. 

Fig. 4. Polarization quantities Do, D n , D p and D q . Lines are the same as in Fig. 3. The 
experimental data in Ref. M are represented by the dots with the error bars. 



